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ABSTRACT 

When a microlensing light curve is contaminated by blended light from unresolved 
stars near the line of sight to the lensed star, the light curve shape and corresponding 
parameterization for the event will differ from the values expected when the event is 
not affected by blending. As a result, blending makes it difficult to identify the major 
lens population and to estimate the amount of lensing matter. In order to estimate the 
effect of blending on the result of lensing experiments, it is, therefore, essential to know 
how the observed lensing parameters change depending on the fraction of blended light. 
Previously, the changed lensing parameters were obtained with a statistical method 
that not only required a large amount of computation time but also was prone to 
uncertainty. In this paper, we derive analytic relations between the lensing parameters 
with and without the effect of blending. By using these relations, we investigate the 
dependence of the observed lensing parameters on the amount of blended light, the 
impact parameter, and the threshold amplification for event detection. 
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1 INTRODUCTION 

Following the proposal of Paczyriski (1986), searches for 
Massive Astronomical Compact Objects (MACHOs) by de- 
tecting light variations of stars caused by gravitational mi- 
crolensing are being carried out by several groups (MACHO, 
Alcock et al. 1997; EROS, Ansari et al. 1996; OGLE, Udalski 
et al. 1997). 

When an isolated source star is gravitationally mi- 
crolensed, its flux is amplified by an amount 
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where u is the lens-source separation in units of the angu- 
lar Einstein ring radius, #e- The lensing parameters io, /3, 
and 4b represent the time of maximum amplification, the 
lens-source impact parameter in units of 9e, and the Ein- 
stein ring radius crossing time (Einstein time scale), respec- 
tively. These parameters are obtained by fitting model light 
curves in equation (1.1) to the observed one. Among these 
parameters, the Einstein time scale is related to the physical 
parameters of the lens by 
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where M is the mass of the lens, v is the lens-source trans- 
verse velocity, D a i, Di s , and D oa are the separations be- 
tween the observer, lens, and source star, and te = D i9e is 
the physical size of the Einstein ring radius. Since the Ein- 
stein time scale results from the combination of the physical 
parameters of the lens, it is difficult to obtain information 
about these parameters of individual lenses. However, one 
can still obtain information about the major population of 
lenses from the distribution of Einstein time scales by mod- 
eling the distribution and motion of lens matter constrained 
from other types of observation (Han & Gould 1996; Gould 
1996; Zhao, Rich, & Spergel 1996; Han & Chang 1998). In 
addition, from the comparison of the experimentally deter- 
mined optical depth, which is proportional to the summa- 
tion of the Einstein time scales of individual events, i.e. 
r oc y^.tB.j; to the theoretically expected value, experi- 
ments toward Magellanic Clouds allow one to constrain the 
MACHO dark matter fraction in the Galactic halo. 

However, the observed light curves of a significant frac- 
tion of lensing events will differ from the light curve in equa- 
tion (1.1). This is because current lensing experiments are 
being conducted toward very dense star fields such as the 
Magellanic Clouds and the Galactic bulge. While searches 
towards these crowded fields result in an increased event 
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rate, it also implies that many of the observed light curves 
are blended by the light from unresolved stars that are not 
being lensed. When an event is affected by blending, the 
observed light curve becomes 
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where Fo is the unblended and unamplified flux of the lensed 
star and B is the amount of blended light. The problem of 
blending is that it is very difficult to detect the presence 
of a blend, and practically impossible to correct for the ef- 
fect by purely photometrical means (Wozniak & Paczyhski 
[1997], see also Figure 2 for example blended light curves 
and the corresponding degenerate unblended light curves). 
Since blending makes the observed parameters of an event 
differ from the values expected when the event is not affected 
by blending, blending leads to misidentification of the ma- 
jor lens population. In addition, the MACHO fraction de- 
termined from the lensing optical depth will be subject to 
great uncertainty. Throughout this paper, we consistently 
use term 'observed lensing parameters' to represent the lens- 
ing parameters determined from a blended light curve not 
knowing it is affected by blending. 

Blending can be classified into several types depending 
on the origin of the blended light. The first type, regular 
blending, occurs when a bright source star registered on a 
template plate is lensed and its flux is affected by the light 
from faint unresolved stars below the detection limit im- 
posed by crowding. Second, lens blending occurs when the 
origin of the blended light is the lens itself (Kamionkowski 
1995; Buchalter & Kamionkowski 1997; Nemiroff 1997; Han 
1998). Third, amplification-bias blending occurs when one of 
several unresolved faint stars in the seeing disk of a reference 
star is lensed, and the flux of the lensed star is associated 
with the flux from other stars in the integrated seeing disk 
(Nemiroff 1994; Han 1997a; Alard 1997). Finally, if a source 
is composed of binaries and only one of the component is 
gravitationally amplified, the light from the other binary 
component simply contribute as blended light: binary-source 
blending (Griest & Hu 1992; Dominik 1998; Han & Jeong 
1998). 

In order to estimate the effects of the individual types 
of blending on the results of lensing experiments, one has to 
know how the observed lensing parameters change depend- 
ing on the fractions of blended light. One of the methods 
to determine the changed lensing parameters is statistically 
fitting a series of unblended light curves to the blended light 
curve of interest (e.g., Wozniak & Paczyhski 1997). In this 
case, determining the parameters requires a large amount of 
computation time for fitting process. Furthermore, for the 
full analysis of the blending effect, one should repeat the 
same procedure for events with all combinations of lensing 
parameters. In addition, there are inevitable statistical un- 
certainties that accompany the use of this method. Di Ste- 
fano & Esin (1995) attempted to develop analytic formal- 
ism for the treatment of blending. However, their so-called 
"blended Einstein ring radius" quantifies only how the max- 
imum allowed impact parameter of lens-source encounter for 
event detection changes as a function of blended light frac- 
tion (see § 2). Since all information about lenses are ob- 
tained from lensing parameters, especially the Einstein time 
scale, simple analytic relations between the lensing param- 



eters with and without blending will be very useful for the 
estimation of blending effect on the result of lensing experi- 
ments. 

In this paper, we derive these relationships and use the 
relations to investigate the dependence of the observed lens- 
ing parameters on the amount of blended light, the impact 
parameter, and the threshold amplification for event detec- 
tion. 



2 CHANGES IN LENSING PARAMETERS 
DUE TO BLENDING 

The shape of a microlensing event light curve is character- 
ized by its height (peak amplification A p ) and width (event 
duration td), which are dependent on the lensing parame- 
ters: f3 for the height (see equation [2.2]) and £e and (5 for 
the width (see equation [2.5]). Because of the change in the 
peak amplification due to blending, the observed impact pa- 
rameter differs from the value determined from unblended 
light curve. Without blending, the impact parameter is re- 
lated to the peak amplification of the unblended light curve 
by 
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When the event is affected by blending, on the other hand, 
the observed peak amplification decreases into 
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As a result, the impact parameter obtained by ignoring 
blending effect becomes dependent on the amount of the 
blended light by 
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Not only the peak amplification, but also the event du- 
ration is affected by blending. As an event is always am- 
plified by more than 1 by definition, we define the event 
duration by half of the time a given event has an amplifi- 
cation greater than some threshold amplification of A t h,o- 
Q For a given threshold amplification for event detection, 
the maximum allowed impact parameter (threshold impact 
parameter) is determined by 
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For example, if the adopted threshold amplification is 
A t h,o = 3/\/5, the threshold impact parameter becomes 
(3th, o = 1. According to this definition of a detectable mi- 
crolensing event, the source star should enter the Einstein 

* To avoid confusion, we note the difference between the event 
duration, t^, and the Einstein time scale, t-^. For an unblended 
event detected with a threshold amplification of A t h q = 3/V5, 
td represent the half of time for the source to stay within the 
Einstein ring while is the time scale for a source to cross the 
Einstein ring radius. The longest source star trajectory within 
the Einstein ring is the Einstein ring diameter. Therefore, the 
two time scales become equal only when the source crosses the 
center of the Einstein ring, but otherwise is always less than 
t E . 
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Figure 1. Upper panel: The lens system geometry of an un- 
blended lensing event. To be detected a source star should enter 
the ring (dotted line) around the lens (L, marked by 'x') with a 
radius /3 t ^ q (threshold impact parameter). If the adopted thresh- 
old amplification is A^o = 3/\/5, ftfi.o = 1; an d the ring is 
equivalent to the Einstein ring. The straight line (with an arrow) 
represents a source star trajectory. The three points on the source 
star trajectory (marked by numbers) represent the source star po- 
sitions at the moments when it enters the ring, approaches most 
closely to the lens, and leaves the ring, respectively. Lower panel: 
The resulting light curve from the lens system geometry in the 
upper panel. The three points on the light curve (marked also 
by numbers) represent the amplifications at the source star posi- 
tions with corresponding numbers in the upper panel. The peak 
amplification (A p ) and the threshold amplification (j4 t j, q) are 
represented by dotted lines. The Einstein time scale (4e) and the 
event duration (td) are marked by thick solid lines with arrows. 



ring. However, since the threshold amplification can be dif- 
ferent values depending on the observational strategy and 
corresponding event selection criteria, we leave the thresh- 
old amplification as a variable. From the definition of the 
threshold impact parameter as /3f h = /3 2 + (td/tE)' 2 , one 
finds that the duration of an unblended event is related to 
the threshold amplification and the unblended lensing pa- 
rameters by 



td(A t h,o, *e, ft) = t-E (ftth,o — /3 2 ) 



2\l/2 



(2.5) 



In Figure 1, to visualize the definitions of ftth,o, A v , and td 
and the relation between tE and td, we present the geometry 
of a lens system for an unblended event (upper panel) and 
the resulting light curve (lower panel). In the upper panel, 
the source star trajectory is represented by a straight line. In 
the lower panel, the peak amplification (A p ) and the thresh- 
old amplification (Ath.o) are represented by dotted lines. 



The Einstein time scales (4e) and the event duration (td) 
are marked by thick solid lines with arrows. 

On the other hand, if the event is affected by blending, 
the same event must have a higher amplification because 
blending increases the threshold amplification to: 



A th (Ath,o, B/F ) = A th , (l + B/Fo) - B/F . 



(2.6) 



Blending also decreases the observed duration of the event 
to: 
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is the decreased threshold impact parameter caused by the 
increased threshold amplification A t h in the presence of 
blending. We note that the formula of so-called "blended 
Einstein ring" of Di Stefano & Esin (1995) is identical to 
our threshold impact parameter in equation (2.8) for a spe- 
cial case for A t h.o = 3/v5. 

As a result of the changes in the impact parameter and 
the duration of event, the Einstein time scale also changes. 
Not knowing an event is affected by blending, one will ob- 
tain the Einstein time scale from the relation in equation 
(2.5), but based on the observed impact parameter and the 
duration of event by 
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In the equation, the threshold impact parameter /3th, o is 
included instead of the decreased value of ft t h due to blend- 
ing, because one still thinks he (or she) consistently applied 
the same threshold impact parameter of /3th, o- Since the ob- 
served duration of the event is related to the unblended lens- 
ing parameters by equation (2.7), one obtains the relation 
between the Einstein time scales with and without blending 
by 



tE.ob 
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ft 2 
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(2.10) 



where the definitions of the various types of impact param- 
eter of ft, ftobs, ftth,o, and ftth are given in equations (2.1), 
(2.3), (2.4), and (2.8). 

In Figure 2, we present several illustrative light curves 
with blending (solid curves) and the corresponding degen- 
erate unblended light curves (dotted curves) whose lensing 
parameters (i.e., /3 f, s and tE.obs) are determined by our de- 
rived relations in equations (2.3) and (2.10). All of these 
events have the same unblended impact parameter /3 = 0.2 
and the adopted threshold amplification is A t h,o = 3/v5 
(represented by a solid line). The fraction of blended light 
and the fractional changes in the lensing parameters with 
and without blending are marked in each panel. As shown 
by Wozniak & Paczyriski (1997), the light curves of blended 
events are very well fit by unblended light curves of different 
lensing parameters. 

With the derived analytic relations, the effects of blend- 
ing on the observed lensing parameters can be computed 
explicitly. We compute these changes in the lensing param- 
eters due to blending with the equations derived above for 
the impact parameter (2.6) and for the Einstein time scale 
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Figure 2. Blended event light curves (solid curves) and the 
corresponding degenerate unblended light curves (dotted curves). 
The lensing parameters determined from the unblended fit to the 
blended light curves (i.e., /3 (, s and tpj f, 3 ) are obtained by using 
the derived equations (2.3) and (2.10). All of these events have 
the same unblended impact parameter ((3 = 0.2) and the Einstein 
time scale. The adopted threshold amplification is A th g = 3/\/5 
(represented by a solid line). The fraction of blended light (B/Fq) 
and the fractional changes in the lensing parameters {j3 obB / f3 and 
*B o6sAe f° r the impact parameter and the Einstein time scale, 
respectively) with and without blending are marked in each panel. 
One finds that the blended light curves are very well fit by un- 
blended light curves of different lensing parameters. 

(2.10). Figure 3 shows how these parameters vary as a func- 
tion of the blended light fraction for different values of the 
impact parameter. To compute the observed lensing param- 
eters, we assume that the threshold amplification for event 
detection is A t h,o = 3/\/5- In the figure, the relation be- 
comes discontinuous above a certain value of B/Fq because 
events affected by more blending than this amount cannot 
be amplified higher than A t h,o and thus cannot be detected. 
As expected, one finds that the observed Einstein time scale 
decreases and the observed impact parameter increases as 
the fraction of blended light increases. We note, however, 
that although these qualitative trends of the changes in the 
observed lensing parameters due to blending are already 
known (Di Stefano & Esin 1995; Alard & Guibert 1997; Han 
1997b; Goldberg 1998) and numerical quantification based 
on statistical methods was performed (Wozniak & Paczyn- 
ski 1997), our derivation is the first to analytically quantify 
these changes. 

Another interesting result from Figure 3 is that the ob- 
served lensing parameters depend in different ways on B/Fq 
for different values of the unblended impact parameter (3. 
This is because even for events with the same Einstein time 
scale and the same fraction of blended light, the peak ampli- 



Figure 3. The fractional changes in the observed lensing param- 
eters (tE.,obs/~t~E and (3 bs/l3 for the Einstein time scale and the 
impact parameter, respectively) as a function of the fraction of 
blended light (B/Fo) for events with different values of the impact 
parameter. We assume that the threshold amplification for event 
detection is At^ g = 3/%/E. In the figure, the relation becomes dis- 
continuous at a certain value of the blended light fraction B/Fq 
because events more strongly affected by this amount of blended 
light cannot be amplified higher than A t ^ q, and would fall below 
the detection threshold. 

fication and the duration of the event increase as the impact 
parameter decreases. The decrease in iE,o6s and the increase 
in f3 bs become more important for events with smaller im- 
pact parameters, implying that blending is more important 
for events with higher amplifications. 

We also find that the observed Einstein time scale de- 
pends differently on the fraction of blended light for differ- 
ent values of the adopted threshold amplification. Figure 4 
shows how tE,obs changes for different values of Ath,o for 
an example event with /3 = 0.2. From the figure, one finds 
that as the threshold amplification increases, the decrease in 
the observed Einstein time scale becomes more important. 
This is because by increasing the threshold amplification 
the observed duration of an event decreases, while the un- 
blended Einstein time scale remains the same regardless of 
the adopted A t h,o- Unlike the Einstein time scale, however, 
the observed impact parameter has the same dependency on 
B/Fq regardless of the adopted values of A t h,o- This is be- 
cause f3 bs is determined solely from the peak amplification, 
which does not depends on the adopted A t h,o- 



3 CONCLUSION 

We have derived analytic relations for the observed mi- 
crolensing parameters with and without the consideration 
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Figure 4. The dependence of the observed lensing parameters 
(*E,obs/*E and f3 bs/P f° r the Einstein time scale and the impact 
parameter, respectively) on the fraction of blended light (B/Fq) 
for different values of the adopted threshold amplification A t h q. 
This model event has an unblended impact parameter of /3 = 0.2. 
As the value of the threshold amplification increases, the decrease 
in the observed Einstein time scale becomes more important. Un- 
like the Einstein time scale, however, the observed impact param- 
eter has the same dependence on B/ Fq regardless of the adopted 
values of A t ^ q. Note that the four lines in the lower panel coin- 
cide. 
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of blending effect. We have used these relations to investi- 
gate the dependence of the observed lensing parameters on 
the fraction of blended light for events with various values 
of impact parameters under various selection criteria. The 
results of this investigation are as follows: 

(i) The observed Einstein time scale decreases and the 
impact parameter increases with an increasing fraction of 
blended light. Although these trends in the observed lensing 
parameters are already known, this derivation is the first to 
analytically quantify this effect. 

(ii) The observed lensing parameters depend in different 
ways on the fraction of blended light for events with differ- 
ent impact parameters. These changes in the observed lens- 
ing parameters are more important for events with smaller 
impact parameters, implying that the effects of blending be- 
comes more important for events with higher amplifications. 

(iii) The observed Einstein time scale depends on the 
threshold amplification for event detection. With increasing 
value of the adopted threshold amplification, the decrease in 
the observed Einstein time scale becomes more important. 
However, the observed impact parameter has the same de- 
pendence regardless of the adopted threshold amplification. 
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